Dynamical topology and statistical properties of spatiotemporal chaos 
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For spatiotemporal chaos described by partial differential equations, there are generally locations 
where the dynamical variable achieves its local extremum or where the time partial derivative of 
the variable vanishes instantaneously To a large extent, the location and movement of these topo- 
logically special points determine the qualitative structure of the disordered states. We analyze 
numerically statistical properties of the topologically special points in one-dimensional spatiotem- 
poral chaos. The probability distribution functions for the number of point, the lifespan, and the 
distance covered during their lifetime are obtained from numerical simulations. Mathematically, we 
establish a probabilistic model to describe the dynamics of these topologically special points. In 
despite of the different definitions in different spatiotemporal chaos, the dynamics of these special 
points can be described in a uniform approach. 



Spatiotemporal chaos is often mediated by cer- 
tain types of dynamical spatial structures or topo- 
logical features such as defects in defect-mediated 
turbulence. The dynamics of such special entities 
characterizes the disordered states and presents 
a simplified description of spatiotemporal chaos. 
For the more general spatiotemporal chaos that 
are not featured by defects and apparent geo- 
metrical structures, it is a challenge to obtain a 
quantitative mathematical characterization that 
records or preserves qualitatively the geometric 
structures of the complex patterns. In this pa- 
per, we propose to use the topologically special 
points such as extremum and critical points in 
disordered states to characterize qualitative prop- 
erties of spatiotemporal chaos. To a large ex- 
tent, the location and movement of these topo- 
logically special points determine the qualitative 
structure of the disordered states. We calculate 
statistical properties of the extremum and critical 
points with the Brusselator model and the com- 
plex Ginzburg-Landau equation. A probabilistic 
model that can regenerate the simulation results 
is proposed. 
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Spatially extended dynamical systems may exhibit 
spatiotemporal chaos (STC) characterized by a fi- 
nite correlation length in both space and time [1]. 
While self-organization in pattern formation systems has 
been intensively investigated in the last decades, high- 
dimensional disordered states are still poorly understood. 
The main challenge in this field has been to reveal the 
onset of STC [2-4 and to establish methods to char- 
acterize such disordered states [H |5l [6]. In many spa- 
tially extended systems, spatiotemporal chaos is charac- 
terized by the presence of defects. For such disordered 
states generally referred to as defect-mediated turbulence 
(DMT), a statistical description of defect dynamics has 
been used as a unifying approach to characterize such 
STCs [THITI. Numerical and theoretical analyses have 
been carried out in the complex Ginzburg-Landau equa- 
tion [7, 8 , the FitzHugh-Nagumo-type system [9 , the 
Willamowski-Rossler model [10 , the model for rotating 
non-Boussinesq convection [11 , and catalytic surface re- 
actions [12]. Experimental studies of defect statistics 
have been also reported from electro-convection in liq- 
uid crystals [I3l [14 , the catalytic surface reaction [15], 
and the Belousov-Zhabotinsky reaction [16]. 

Similar to the defects in DMTs, other spatial struc- 
tures or topological features have been used to quantify 
the dynamical complex patterns. Such examples have 
been reported from algebraic topology for measuring the 
complexities in disordered states [18 , instantaneous stag- 
nation points with vanishing velocity in the velocity field 
of fluids Ll9j, and self-replicating spots in the Gray-Scott 
model [201 EI]. 

In order to understand and quantify the more gen- 
eral STCs that are not featured by defects and appar- 
ent geometrical structures, we here propose a quantita- 
tive mathematical characterization that records or pre- 
serves the dynamical structures of disordered states. For 
STCs described by partial differential equations, there 
are generally locations ro in the scalar field X(r, t) where 
the dynamical quantity achieves its local extremum, i.e.. 
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^'^dr'^^ |ro = 0, as depicted in Fig. la. These instan- 
taneous maximum or minimum points preserve the ge- 
ometrical structures of the disordered states (Fig. Ic). 
If the locations and movements of all of these special 
points are known, the pattern can be primarily deter- 
mined. Similarly, as in the velocity field of fiuid fiows 
[19], the spatiotemporal chaotic field X(r, t) can be char- 
acterized by locations ro where the time partial derivative 
of the dynamical quantity vanishes instantaneously, i.e., 
^^m*^ |ro = 0, as demonstrated in Fig. lb. These criti- 
cal points carry the bulk of the information contained in 
the disordered patterns (Fig. Id). Such defined topolog- 
ically special points, which are present even in ordered 
states, are distinct from the topological defects in previ- 
ously studied DMTs. They are created and annihilated 
irregularly in pairs, changing simultaneously the topolog- 
ical structures of the disordered states. To a large extent, 
the location and movement of these points determine the 
qualitative structure of the disordered states. 

We analyze numerically statistical properties of the 
topologically special points in one-dimensional STCs 
with the Brusselator model and the complex Ginzburg- 
Landau equation (CGLE). The probability distribution 
functions (PDF) for the number of extremum points and 
that of critical points are obtained from numerical sim- 
ulations. PDFs for the lifespan as well as the distance 
covered during the lifetime are also calculated. Mathe- 
matically, we establish a probabilistic model to describe 
the dynamics of these topologically special points. The 
analysis gives results that are in good agreement with 
numerical simulations. In despite of the different defi- 
nitions in different STCs, the dynamics of these special 
points can be described uniformly. 

Models and numerical results. — The reaction-diffusion 
system with Brusselator reaction kinetics is a standard 
model system for the investigation of self-organized struc- 
tures in nonlinear chemical reactions, which is governed 
by the following set of partial differential equations [22^, 



dtX = Dxd^X 
dtY = DydlY- 



X^Y- 



-(^ + 
BX. 



1)X + A, 



(1) 
(2) 



Figure 1 demonstrates the chemical turbulence generated 
by Eqs. 1 and 2 with parameters A = 2.0, B = 5.5, Dx = 
1, Dy = 0. Figure la and lb depict the snapshots of the 
one-dimensional fields X(r, t) and ^^^p^, respectively. 
The maximum (minimum) points in Fig. la are denoted 
as red (blue). The dots in Fig. lb are for the critical 
points where ^^qI'^'' = 0, with the red or blue for loca- 
tions where dr[dtX{r,t)] > or dr[dtX{r,t)] < 0. As 
these special points demonstrated in Fig. la and lb are 
always created and annihilated in red-blue pairs, the red 
and the blue can be denoted with positive (+) and nega- 
tive (-) signs similar to the positive and negative charges 
for topological defects in DMTs. The space-time plots 
in Fig. Ic and Id demonstrate the time evolutions of 
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FIG. 1: One-dimensional spatiotemporal chaos generated 
from the reaction-diffusion system with Brusselator reaction 
kinetics (Eqs. 1 and 2). (a) Snapshot of X(r,t) and local 
maximum (red) and local minimum (blue) points; (b) Snap- 
shot of dtX(r, t) and critical points where dtX{r^ t) = 0, with 
red for dr[dtX{r,t)] > and blue for dr[dtX{r,t)] < 0. (c) 
and (d) are space-time plots (horizontal space and vertical 
time) depicting the time evolutions of the topological special 
points in (a) and (b), respectively, (e) The space-time map 
of X(r, t) coded in gray scales. The arrows in (a) indicate the 
red-blue pairs being created or annihilated in the disordered 
states. Parameters: A — 2.0, B — 5.5, Dx — l^Dy — 0. 
Periodic boundary conditions are used in simulations. 



the topologically special points in Fig. la and Fig. lb, 
respectively. The gray-scaled space-time pattern in Fig. 
le is generated from X(r, t). To a large extent, the com- 
plexities of the disordered states (Fig. le) are primar- 
ily preserved or recorded in the time evolutions of these 
topological points (Fig. Ic and Id). The definitions of the 
special points for the Brusselator can be easily extended 
to the CGLE and other systems. The CGLE which is 
one of the most-studied nonlinear equations in physics is 
given by [23 , 

dtA = A + (1 + ici)d^A - (1 - ics)\A\^A. (3) 

Beyond the Benjamin- Feir instability, three types of 
chaotic behavior are possible in the one-dimensional 
CGLE, namely phase turbulence, space-time defect tur- 
bulence [24^, and intermittency [25|. As the emphasis of 
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FIG. 2: (a) Time series N(t) for the number of extremum 
points, (b) The normaUzed pair correlation function a(r) for 
the extremum points (inset for the critical points), (c) Bi- 
modal probability density function p{x,t = 2.5) for extremum 
points (inset for critical points). The arrows indicate the di- 
rections p{x, t) is splitting. All the simulation results are ob- 
tained from the Brusselator model with the same parameters 
in Fig. 1. 



this paper is on the presentation of the technique, we here 
concentrate on the topological special points as defined 
for the Brusselator without discrimination in the types 
of STCs. 

As a simplified representation of spatiotemporal chaos, 
the dynamical behavior of extremum and critical points 
are more convenient to describe and cope with. The 
points perform diffusive and drift motions in the space 
after they are created in pair, are destroyed when two 
neighboring points with opposite sign bump into each 
other. The processes occur continuously and indefi- 
nitely, and the spatiotemporal pattern looks like a chaotic 
"soup" . The number of extremum or critical points in the 
disordered pattern undergoes fast fluctuations. Figure 2a 
illustrates a time series N{t) of the number of extremum 
point calculated from the Brusselator. The fast stochas- 
tic fluctuations in N{t) is a manifestation of the disorder 
and incoherence in the turbulence. Figure 2b demon- 
strates the pair correlation functions for the extremum 
and critical points, showing short range interactions be- 
tween these points. 

The diffusive and drift motion of extremum and critical 
points can be examined from the time dependent prob- 



FIG. 3: Numerically simulated gain and loss rates, which are 
approximately linear functions, for extremum points obtained 
from the Brusselator model (a) and critical points from the 
CGLE (b), and probability distribution functions for the num- 
ber of extremum points obtained from the Brusselator (c) and 
of critical points from the CGLE. Theoretical PDFs from Eq. 
4 are depicted as solid lines. 



ability density function p(x, t) for the distance x that a 
point travels from the location where it is created at the 
time t from the moment it is created. This can be cal- 
culated by tracking the points from their creation to an- 
nihilation processes. Apparently the initial distribution 
p(x,0) is a delta function S{x). As demonstrated in Fig. 
2c, it will then split as the pairs of point that are newly 
born with opposite signs move apart, forming a bimodal 
distribution p{x^t). By monitoring the creation and an- 
nihilation events, the gain and loss rates can be obtained 
numerically (Fig. 3a and 3b). Both are approximately 
linear functions in the number of points in the system. 
From the time series N(t), the probability density func- 
tion (PDF) P{N) which is the probability for finding N 
points in the disordered states can be derived. Figure 3c 
and 3d (dots) demonstrate the PDFs for the extremum 
and critical points obtained from the Brusselator and the 
CGLE, respectively. 

By monitoring the creation, movement, and annihila- 
tion processes of these points, we derive two extra prop- 
erties: the probability density function P(t) for the life 
span T of the special points, and the distribution P{X) 
of the displacement A (i.e., the distance a point covers 
during its lifetime). Figure 4 (dots) demonstrates the 
life span distributions for extremum and critical points 
calculated separately from the Brusselator model (Fig. 
4a) and the CGLE (Fig. 4b). They are uniformly expo- 
nential decays. Figure 5 (dots) depicts the distribution 
of displacement A for extremum and critical points in 
the Brusselator model (Fig. 5a) and CGLE (Fig. 5b), 
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FIG. 4: Lifetime distribution P{r) of the topologically special points obtained from numerical simulations (black dots) and 
from theoretical prediction of Eq. 9 (red line), (a) Extremum points in the Brusselator model, C = 0.7482; (b) Critical points 
in the CGLE, C = 0.1964. (c) Critical points in the autocatalator model, C = 0.09669. 



respectively. The distributions are also uniformly tent- 
shaped. The left-right symmetry in the distributions is 
due to the symmetrical splitting movement of the special 
points after they are created. 

Theoretical analysis. — We now seek to establish a 
probabilistic model to generate quantitatively the numer- 
ical results. The topologically special points we defined 
here are similar to defects in DMTs. The master equa- 
tion description of defect dynamics [7] can be extended 
in straightforward to the situation we consider here. The 
probability distribution function for the number of ex- 
tremum or critical point can be readily expressed as [7], 



P{N) 



p{o) n 

j=0 



(4) 



where 11+ and are, respectively, gain and loss rates 
that define the transition rates of the Markov chain. 
From our simulations, S+ and I]_ for the Brusselator 
and the CGLE are well approximated by linear depen- 
dence on AT, S+(Ar) = co+ci7V, S-(7V) = ao^aiN (solid 
lines in Fig. 3a and 3b). The resulting distributions from 
Eq. 4 are depicted in Fig. 3c and 3d (curves) which are 
in good agreement with simulations in despite that the 
special points are different and in different models. 

Our finding that the creation and annihilation rate of 
the special points are both linear with N is distinct from 
previous findings in 2D DMT where the annihilation rate 
is proportional to A^^ [7j. The difference might be due 
to that different dimensions have been considered. The 
N'^ annihilation rate is intuitive in 2D systems. In ID 
systems as we consider, the positive spots and negative 
spots are alternatively distributed in the ID space. As 
an elementary annihilation or creation event occurs be- 
tween two spots with opposite signs, it is correlated with 
two neighboring spots. The annihilation/creation rates 
should thus be proportional to the number of spot pairs 
in the space and would result in the linear dependence 
on N as we observed. 

In order to analyze the behaviors of topologically spe- 
cial point in more detail, an equation for the probabil- 
ity density function p(x, t) needs to be established. As 
observed in our simulations, the extremum and critical 



points perform drifting and diffusive motion (Fig. 2c), 
and suffer random extinctions during these processes. To 
describe these properties, we propose a modified Fokker- 
Planck equation for the time evolution of p+(x,t) and 
p-{x,t) with p{x,t) = {p+{x,t) -|-p_(x,t))/2, 

X — Xi{t))]}p±{x,t). 

i 

(5) 

The sign ± describes the symmetrical splitting of p(x, t) 
into p+ and p_. In Eq. 5, constant drifting and dif- 
fusion are assumed. The right most term is used to 
model interaction decays, where Xi{t) is the location of 
a special point leading to the decay. In a mean-field ap- 
proximation, the last term is proportional to the nor- 
malized pair correlation function a{x) (Fig. 2b). One 
has < c^-5{x — Xi{t)) >t= c(j{x). Due to that (j{x) 
approaches a constant value beyond a small spatial cor- 
relation in STC, one has. 



The above equation has the solution 

p±ix,t) 



1 (xTvtf. . . 
^^Pi 77^::;^ — J exp(— Ct). 



4Dt 

The probability p(x, t) is thus, 

p(x,t) = (p+(x,t)+p_(x,t))/2. 



(6) 



(7) 



(8) 



The life span distribution P{r) can be calculated directly 
from the decreasing of the total probability. 



P(t) 



_d_ 
dr 



f 

J — c 



p{x,r)dx = Cexp(-Cr). (9) 



The probability is exactly an exponentially decaying 
function where C depends on the detailed short-range 
interactions. This theoretical result is in good agreement 
with numerical results as shown in Fig. 4a and 4b. 

To calculate the distribution of displacement P(A), we 
have the probability density for the annihilation at loca- 
tion X with life time r as pF{x^r)P{r)^ where ppix^^t) — 
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FIG. 5: Displacement distribution P(A) for the special points obtained from numerical simulations (black dots) and from 
theoretical predictions of Eq. 11 (red line), (a) for extremum points in the Brusselator model, v = 4.331, = 2.9925; (b) for 
critical points in the CGLE, v — 0.764, D — 0.0601; (c) for critical points in the autocatalator model, v — 3.0, D — 2.526. The 
values of C are as in Fig. 4. 



(p+(x,t) + p_(x, t))/2|c^o is the solution of Eq. 6 with 
C = 0. The free path length distribution P(A) is, 

nOO 

P(A)= / dTpF{\T)P{T). (10) 

The integration leads to, 

= 7ifc^^P[-2^(-^^^C^T^)^] 

x(exp[-vA/2D] +exp[vA/2D])/2. ^ ^ 

For a specific disordered state, the drift velocity v and 
diffusion coefficient D in the above equation can be de- 
termined from the simulations of p(x, t) (refer to Fig. 2c) 
using Eq. 8, and the constant C can be determined from 
the simulated life span distribution according to Eq. 9. 
With these parameters, the theoretical result predicted 
by Eq. 11 for the extremum and critical points the Brus- 
selator and in the CGLE demonstrated in Fig. 5a and 
5b (curves), respectively, which also agree uniformly well 
with the results obtained from direct simulations. 

Summary and Conclusion. — In summary, we have 
proposed using topologically special points such as ex- 
tremum and critical points in disordered states to charac- 
terize the qualitative properties of spatiotemporal chaos 
that is not mediated by defects. The special points can 
be denoted as positive or negative, and undergo pairwise 
interactions, changing the topological structures of the 
disordered states when they are created and annihilated. 
The dynamics of such points are similar to the A-\-B ^ C 
type chemical reaction in ID space ^26|, where different 
reactants A and B combine into resultant C when they 
encounter and C decomposes into A and B. We have 
checked the dependence of the average on the length of 
the system (data not presented here). The average num- 
ber of special points is proportional to the size of the sys- 
tem and is therefore extensive in spatiotemporal chaos. 

We demonstrated with the Brusselator model and the 
CGLE that the dynamics of special points determining 
the qualitative properties of STCs can be described in a 
probabilistic model, in spite of different definitions of the 



special points and different dynamical systems. Specifi- 
cally, the previous statistical description of defect dynam- 
ics in defect-mediated turbulence was extended to the 
general spatiotemporal chaos. A modified Fokker-Planck 
equation was proposed as a uniform description for the 
drifting, diffusive and extinctions processes. Theoretical 
analysis gives explicitly expressions for the lifespan dis- 
tribution and the displacement distribution, and predicts 
correctly numerical simulation results. 

The definition of the topological and dynamical spe- 
cial points and the probabilistic approach that we pro- 
posed for ID STCs in the Brusselator and CGLE can 
be extended to other disordered states. For instance, 
we have analyzed the chemical turbulence which emerges 
from the instability of a traveling wave. The model is 
a two-variable, cubic autocatalator with equal diffusiv- 
ities of the species which is described by the following 
equations [ 27] , 

dta = 5dla^l- a- iiap'^, (12) 
dtp = dlp^iiap'^ -(1)13. (13) 

As demonstrated in Fig. 4c, the lifetime distribution of 
the critical points obtained from numerical simulations 
(black dots) in the STC regime, with (5 = 1,^ = 3.0, ji = 
38.0, is in satisfactory agreement with that predicted by 
Eq. 9 (red line). Figure 5c depicts the same agreement 
between the numerical simulation and theoretical predic- 
tion of the displacement distribution. 

In 2D defect-mediated turbulence, the defect is a con- 
venient characterization of turbulence as previously in- 
vestigated. The method present here is suitable only for 
ID systems. It is convenient to define and handle topo- 
logically special points in ID systems. In 2D systems, 
the special points as we defined in ID would be no longer 
isolated points, but could be connected into ID lines or 
filaments. 

The work is financially supported by the NSFC 
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